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Abstract 



■ , We derive the effective Lagrangian for the low-energy massive meson excitations 

D , of the color-flavor-locking (CFL) phase of QCD with 3 flavors of light quarks. We 

compute the decay constants, the maximum velocities, and the masses of the mesons 
at large baryon chemical potential ^. The decay constants are linear in /x. The meson 
maximum velocities are close to that of sound. The meson masses in the CFL phase 
are significantly smaller than in the normal QCD vacuum and depend only on bare 
quark masses. The order of the meson masses is, to some extent, reversed compared 
to that in the QCD vacuum. In particular, the lightest particle is r/'. 



1 Introduction 

The behavior of QCD at finite baryon density could affect the physics of neutron stars, 
supernovas and of heavy-ion colhsions. Quark pairing at the Fermi surface leading to diquark 
condensation and color superconductivity is a subject of many recent theoretical studies pl- 



15| . In particular, as first pointed out by Alford, Rajagopal and Wilczek g], in the case of 
3 light fiavors the diquark condensate can "lock" color and fiavor symmetry rotations. The 
result is the color-fiavor-locking (CFL) phase where an interesting pattern of chiral symmetry 
breaking emerges: SU{3)l x SU{3)r x S'f/(3)coior -^ SU (3) l+r+coIot- The CFL phase has 
many similarities with the chiral symmetry breaking QCD vacuum and nuclear matter. This 
observation has lead to the conjecture that quark and nuclear matter might be continuously 
connected 0. 

As already pointed out in Refs. 0, H, the CFL phase is characterized by a hierarchy of 
energy scales. At the lowest scale lie 10 (pseudo-)Goldstone modes, arising from the breaking 
of axial fiavor symmetry, baryon U{1)b symmetry, and axial U{1)a symmetry. These modes, 
except for the ^7(1)^, would be massless if quark masses were zero, but in reality they do 
have small masses. At a higher scale there are quark excitations, which are separated by 
the superconducting (BCS) gap A, and the gluons which acquire mass of order gfi from 
the Meissner effect. The Goldstone modes, therefore, dominate the physics at energy scales 
smaller than A and can be described by an effective theory. When quarks are massless, the 
Lagrangian of such a theory can be shown to be just the QCD chiral Lagrangian, which 
contains the decay constants and the velocities of the Goldstone bosons as free parameters 
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In this paper, we show that all parameters of the chiral Lagrangian, including the mass 
term, can be completely determined in the weak-coupling regime, or the regime of very high 
densities. Our results can be summarized as follows. We will use the "vacuum" notations 
for the Goldstone bosons, where we have an octet of pseudoscalars vr, K and r] and a singlet 
pseudoscalar rj'. The Goldstone boson related to the breaking of U{1)b will be denoted as 
H. To the leading order in strong coupling, the decay constants of particles in the octet and 
of the singlets i]' and H are linear in the chemical potential, /i: 

^^0.209/i and /^ = /^ ~ 0.195/i. (1) 

The dispersion relation of the mesons has the form e^ = v'^p'^ + m'^, where the maximum 
velocities of all the mesons are close to the speed of sound in ultrarelativistic fiuids: 

1 

V-k = Vh = Vrj' = -l=. (2) 

The H meson remains exactly massless, while all the others acquire masses when nonzero 
quark masses rriq are taken into account. Due to the approximate U{1)a symmetry at large 



fx the meson masses are not proportional to my^, but to rrig 0. In the hmit muim^ -C rris, 
5 mesons (r/ and the kaons) have masses of the order m^, while 4 others (r/' and the pions) 
have masses of the order {mu^dXrisY''^. The masses of the charged pions and the kaons are 
given by: 

ml± = C{mu + md)ms + 2C"(m„ + md)(2m„ + 2md + rris); 
m|-o = C{ma + ms)mu + 2C'{md + ms){2ma + 2ms + "^«); 
mj^± = C{mu + ms)md + 2C"(m„ + m^) (2m„ + 2m^ + m^) ; (3) 

where C ~ 1.578, and C ~ 0.04216. The states of 7r°, r] and 77' mix and their mass matrix is 
complicated. In the limit tUs ^ fnu^d the heaviest is the r] meson, while the rj' is the lightest. 
As we shall see, the meson spectrum of the CFL phase, although similar to that of QCD 
vacuum, has a certain inverse mass ordering: in particular, the neutral kaons are lighter 
than the charged kaons, and the lightest particle is the rj' singlet. The kaons would also be 
lighter than the pions if it was not for a seemingly small, but important, second term in (|^), 
proportional to C . We shall give a simple explanation of these facts below. 

2 The effective Lagrangian 

In this section we review the basic arguments of Ref. [0. The ground state of the CFL 
phase is characterized by the following diquark condensates^ 

^ia _ g*jfc^afec^^6i^cfc^* and F*" ~ e^^V"=(V^^^'V'R )*, (4) 

where the complex conjugation was added for convenience so that X and Y transform under 
SU{'i)c X Sf/(3)L X SU{'i)R as (3,3,1) and (3,1,3) respectively: 

X -^ UlXUJ and Y -^ UrYUJ. (5) 

The low-energy excitations in the CFL phase are given by the slow rotations of the phases 
of X and Y. Therefore, we can factor out the norm of the condensates and consider unitary 
matrices X and Y. Together they give us 9 + 9 = 18 degrees of freedom. 8 of them are 
eaten by the gluons through the Higgs mechanism, and the surviving 10 become low-energy 
excitations. 

One of the Goldstones corresponds to the spontaneous breaking of the U{1)a symmetry, 
which, being a symmetry of the Lagrangian, is violated at the quantum level by the axial 
anomaly, or, equivalently, by instanton-induced interactions. This violation, however, be- 
comes very small at large fi due to the effect of screening which suppresses the instanton 

^There is also an admixture of 6-plet condensates, which, however, is small M. 



density by a high power of l//i [|T^. Therefore, we can treat the U{1)a on equal footing with 
other global symmetries in this regime. 

For simplicity, we shall at first ignore the overall U{1) phases of X and Y and assume, as 
in Ref. |]13|, detX = detF = 1. The Lagrangian should be symmetric under the SU{3)c x 
SU{3)l X SU{3)r rotations @. This condition fixes the Lagrangian to the leading (second) 
order in derivatives: 



C 



eflf 



^Tr 
2 



{X^doXy + [Y^doYY] + spatial gradients 



(6) 



The spatial gradients enter in a similar way but with a different constant instead of /vr. The 
cross term Tt{X'^ d^X) (Y^^ doY) |T3| is allowed by the symmetries, but, as we shall see in Sec. 
H, it is suppressed in weak coupling (i.e., at large /i.) Roughly speaking, to leading order in 
g"^, left and right quarks decouple from each other. 

Since SU{3)c is a local gauge symmetry we must replace the derivatives in Eq. (|) by 
covariant derivatives, 

J 77 



c 



eflF 



-Tr 

J TT 



+ ... 



{XdoX^ - gAof + {YdoY^ - gAoY 

{XdoX^ - YdoY^f + {XdoX^ + YdoY^ - 2gA,f 



Tr 



+ 



(7) 



The second term is responsible for the Higgs effect: the vector-like fluctuations, dX = dY, 
of SU(3) phases become longitudinal components of the gluon A^. The gluons acquire a mass 
of order 0{gf.n) (electric and magnetic masses are, in general, different). We shall see that 
/vr ~ /W, and thus the gluon mass is much larger than the momentum scales p that we are 
considering (p < 2 A -C g^), so gluons decouple from the low-energy theory. The axial-like 
fluctuations of the phases, dX = —dY, can be written as fluctuations of the phases of a new 
unitary matrix: 

S = XY\ (8) 



and the effective Lagrangian takes the form (in Euclidean space): 

J TV ' 



c 



eflF 



-Tr 



daEdnE^ + Vldim^^ 



(9) 



which is the usual Lagrangian of the nonlinear sigma model except that the speed of the 
mesons, f^r, can be different from the speed of light. The matrix S is a singlet in color and 
transforms under SU{3)l x SU{3)ji as 



and describes the meson octet. 



(10) 



What happens if we take into account the U{1) phases of X and Y7 Then it is possible 
to add into the chiral Lagrangian a term proportional to 

{TiX^doXy + (JiY^doYy. (11) 

It is possible to add a cross term, (TTX^doX){TTY^doY) but it is suppressed in weak coupling 
by the same reason by which the cross term was omitted in Eq. (P). We shall make the 
consequences of the term (^) clear by splitting X and Y into SU{3) and U{1) parts, 

X = Xe^'^+^'^ and y = Ye-^'^+^"^; (12) 

where X, Y are SU{3) matrices, and the angle is the variable conjugate to the baryon 
charge, normalized to 1 for a single quark. The normalization of the f/yi(l) phase 6 is fixed 
analogously. Consequently, the field S defined in Eq. (^) now has the form 

S = S e^'^, E = XY^ . (13) 

In terms of S, and 6, the lowest order chiral Lagrangian, consistent with the symmetries, 
has the form 

f 



^Ceff = ^Tr 



d,td,t^ + vldMt^ +12/3 {d,ef + vl,{d,e) 



+l2ffj[{doct>Y + v'H{d4Y\- (14) 

The difference between the decay constants of rj' and H mesons from /,r arises from the 
additional allowed term (p!T|), while the difference between f^i and fn arises from the cross 
term (TrX^9oX)(Try'"9oF), and, therefore, is small at weak coupling. We shall use fr,' = fn 
in the rest of the paper. 

The elementary meson fields vr"^, r]', and H, are defined as 

S = exp(2— ^), e = , 0=—=—, (15) 

\ U J V24U y^fn 

where A"^ (A = 1 ... 8) are Gell-Mann matrices normalized so that TrA'^A'^ = 25^^ . 



3 Decay constants of the Goldstone bosons 

Let us now show that all the decay constants of the Goldstone bosons, /^r, /,,', and Jh-, can be 
computed in the high-density regime and are all proportional to the chemical potential. We 
will demonstrate our method on fu-, whose calculation is the simplest. Let us imagine that 
baryon symmetry is not a global symmetry, but instead a local one. So we introduce into 
the theory a gauge field A^ coupled to the baryon current. Due to the Higgs mechanism. 



the gauge bosons acquire a finite mass, which can be computed in both the microscopic 
theory, where it is expressed via the chemical potential, and in the effective theory, where it 
is proportional to fn- Comparison of the two results leads to the determination of fn- 

In order to compute fn we have to deal only with the part of the effective chiral La- 
grangian ( p^lf ) that contains the phase 0. Since (p is not neutral with respect to the baryon 
symmetry, which is now local, one should replace the derivatives by covariant ones, 

C = I2fl [(9o0 + eA,f + vl{d4 + eA,f] , (16) 

where e is some arbitrary small coupling constant. From Eq. (0) we find that the gauge 
bosons acquire finite masses, which are different for Aq and Aj, 

m% = 2Ae^fl (17) 

m% = 2Avle'fl (18) 

On the other hand, the masses in Eqs. (|T^,|T^) can be computed from the microscopic 
theory, where m^^ has the meaning of the Debye mass (or the inverse Thomas- Fermi screen- 
ing length) of the electric A^ field, while m\_ is the Meissner mass (or the inverse London 
penetration depth) of the magnetic components of A^. To compute these masses, it is most 
convenient to use a low-energy effective theory containing only fermion modes near the Fermi 
surface, similar to the one derived by Hong [0]. We shall work in Euclidean space where 
the action has the form 

^ = /^(^^(p)(^Po + 6p)^(p) + e^t(p)^(p)(,^^(0) + ^^A,(0)))+|d^x^A,A,(19) 

where Vi = Pi/\p\, and we have kept only the gluon modes with zero momentum in the 
interaction term since only they are needed in future discussion. By using the effective 
Lagrangian (|T^ one can avoid dealing with the Dirac structure of the quark propagator, 
which is somewhat complicated in the superconducting phase. The "bare Meissner mass" 
term proportional to rriQ in Eq. ([19|) emerges from integration out of all degrees of freedom 
except the ones near the Fermi surface; its magnitude can be found from the following 
simple argument. In the first-quantization picture, the Hamiltonian of a relativistic particle 
is simply H = |p|. As we couple the particle to the gauge field A^, its Hamiltonian becomes 

e2 , , 

Ti: = |p + eA| + leAo = |p| + ieAo + e{v ■ A) + -— (A^ - (v • A)^ + . . . (20) 

z|p| ^ ^ 



In the second quantization language, the first three terms in Eq. ( pof ) correspond to the first 



term of the Lagrangian (|19|). The last term, if we sum over all particles with momentum p 





Figure 1: The leading order diagrams contributing to decay constants. 
inside the Fermi sphere, reproduce the bare Meissner mass term in the effective Lagrangian 



m^ = 6e^A_ (21) 



(|T|) with 

2 

A consistency check for the Lagrangian (|19]) is the computation of the physical Meissner mass 
in the normal phase. It can be shown that the bare value (^) is exactly the one required 
for the Meissner effect to be absent in the normal phase. 

In the CFL phase, the quark propagator has its simplest form in the basis 

^a^ = E ^r (22) 

A=l V 2 



where A"^ are Gell-Mann matrices ii A = 1 ... 8 and A^ = a/2/3. The Nambu-Gorkov quark 
propagators are diagonal in this basis 

SAB^A 

plT^ + Ay 
{r{pW{-P)) = \ ,"7 ,'7/ ; (24) 



^1 
XAB/\/ 

{^^{P)r{-P)) = „, , ,, , ,, ; (23) 

p 



pl + 4 + Ay 



where e^ = |p| — yU is the energy relative to the Fermi surface. The gaps A^ with A = 1 . . .8 
are equal, but different from Ag. In weak coupling, if we denote Ai = ■ ■ ■ = Ag = A, then 
Ag = -2A.g 

To the leading order, the only contribution to Tn\g is from the sum of two one-loop 
diagrams as shown in Fig. (|ID, which, for zero external momentum, is equal to 

^fi^ f dpodep 9 . (ipo + ep)^ A^ \ 

2 J (27r) tS iPl + ^l + A\r iPl + 4 + ^Ayr ^ ' 

where we use the formula for the phase space near the Fermi surface, / 7|^ = ^ /^ dep. 
The overall factor of 2 in Eq. (^5]) comes from summation over left- and right-handed quarks 

^Our A"* with A =^ 9 and A^ corresponds to Ag and Ai in the notation of Ref. [Q. 



in the internal loop. While the second diagram in Fig. (|l|) is finite, the first, treated formally, 
has logarithmic divergence. The prescription to deal with the divergence is known |16| , and 
is essentially to perform the integration over pq first. We find that the two diagrams in Fig. 
(0) give equal contribution to the squared Debye mass which is given by 

ml = 18e^ 1^. (26) 

This value is the same as the Debye mass that A^ would have in the absence of the su- 
perconductivity. The origin of this coincidence can be made clear by the following argument. 
What we have computed is simply the 00 component of the polarization operator of A^, IIoo- 
Since Aq is coupled to the baryon charge n = ip^tp, Hqq has the following interpretation in 
the linear response theory: the baryon charge density generated by an external uniform field 
Aq is given by 

n = noo(0)Ao. (27) 

But the uniform Aq field is simply a shift of the chemical potential, or the Fermi energy. 
Therefore, IIoo is equal to e^ dn/dfi, i.e., the density of states near the Fermi surface, which 
is exactly the right hand side of Eq. (|26|) . It is easy to see that dn/dfi is the same in the 
normal phase and superconducting phase, provided that the gap in the latter is small. 

Comparing Eq. (p6|) with what one gets from the effective theory, Eq. (|l3), one finds the 
decay constant fn, 



4 27r2 



fn = -.^- (28) 



An important remark is in order here. Eq. (^) tells us that fn depends only on the chemical 
potential yU, but not on the gap A. This might seem contradicting the fact that at A = the 
U{1)b symmetry would be restored and no Goldstone boson is expected in the theory. The 
explanation is that, as A decreases, the domain of applicability of the effective Lagrangian 
p < 2 A shrinks and disappears at A = 0; therefore the persistence of fu does not contradict 
the restoration of U{1)b symmetry at A = 0. 

To find the Meissner mass, we have to evaluate the same two diagrams of Fig. (|1]) , where 
the vertices are attached to Ai and Aj. Instead of being equal, the two diagrams now cancel 
each other, since the vertex factor of the first diagram is ViVj, while in the second diagram 
it is Vi{—Vj). Therefore, the Meissner mass in the superconducting phase is equal to the 
bare Meissner mass in the effective Lagrangian ([191) . Compared to Eq. ([ISD and taking into 
account Eq. (H), one finds that the velocity of the H boson is equal to l/\/3, i.e. the speed 
of sound in relativistic fiuids, 

vl = \. (29) 




Figure 2: A higher-order diagram contributing to the decay constants. One loop contains 
left-handed quark, the other has right-handed quark. 

This result is not completely surprising, since the U{1)b phase (p of the condensate is the 
variable conjugate to the baryon density ip'^ip, whose fluctuations give rise to the sound. 
Therefore, H quanta can be considered as phonons. But they are not hydrodynamic phonons 
since they exist outside the hydrodynamic regime (as in our case, at zero temperature where 
the mean free path diverges.) Moreover, as we shall explain, if one increases the temperature 
the velocity of the Goldstone bosons decreases and becomes zero at the critical temperature, 
while the speed of the hydrodynamic sound is almost insensitive to the temperature in this 
range. Therefore, the fact that the velocity of H quanta (and other Goldstone bosons, as 
we shall see) at zero temperature is equal to the speed of hydrodynamic sound should be 
considered as a coincidence bearing no fundamental reason. 

For the 77' meson, one can repeat the same calculation and see that, to the leading order, 
all diagrams remain the same. Thus we find that /^/ = fn, and that the rj' meson also 
propagates with the sound speed. There is, however, no symmetry that requires the decay 
constants of H and 1]' to be equal; in fact, they are not equal in higher orders of perturbation 
theory. For example, the digram drawn in Fig. ^ gives contribution of opposite signs to fn 
and /,,/. 

The computation of /^ and v,, is completely analogous to that of fn- We could introduce a 
fictitious field coupled to the fiavor currents, but we can also make use of the existing couphng 
to the gluons to compute /^r. The only additional complication is that the interaction vertex 
now has a non-trivial structure. In our basis (^2]) , the interaction vertex has the form 



2 E E tr(^^^''^'')^^^^''^'' (30) 

^ A,C=IB=1 

Straightforward calculations show that the first diagram in Fig. |l| gives (3/4)5f^/i^/(27r^) and 
the second diagram contributes — ((3 + 41n2)/18)(7^/i^/(27r^) to the Debye mass. Taking into 



account the factor of 2 from left- and right-handed quarks, the result reads 

2 21 -8 In 2 (7V ,^^. 

rriA = — (31) 

^° 18 27r2 ^ ' 

The mass in Eq. ( PT| ) is not equal to the Debye mass in the normal phase, in contrast to the 
previous case of the fictitious U{1) boson coupled to the baryon current. 

Notice that if the gluons are coupled to the left-handed quarks alone, the squared Debye 
mass would be half smaller, which is exactly what one obtains by throwing away Y from 
the Lagrangian (0). However, if we add the XY cross term to the Lagrangian, this will no 
longer be true. Therefore, this cross term, though not forbidden by the symmetry, has a 
small coefficient in weak coupling. This can be seen also from the fact that the coupling of 
the left- and right-handed flavor currents is zero at the leading order. In higher orders of 
perturbation theory it receives contribution from the diagram like in Fig.|], where one vertex 
corresponds to the left-handed flavor current and the other to the right-handed one. 

In the effective theory, the gluon mass is given by m\^ = g'^f^. Therefore, one can 
determine the decay constants of the mesons in the pseudoscalar octet. 



_ tlZJll^^ (32) 



21-81n2 /i2 

Ts 2^2 

The ratio f^/f^, = fl/ fl = 2(21 - 81n2)/27 ^ 1.14 is quite close to one. This fact seems 
to stem from the OZI rule |]T^ . 
The Meissner mass is equal to 

^" 27r2V2 4 54 / 54 27r2 ' ^ ' 

where, in the parenthesis in the intermediate expression, the flrst term comes from the bare 
Meissner mass, while the two last terms come from the two diagrams in Fig. 0. Again, in 
contrast with the case of U{1)b, the contribution of the two diagrams does not cancel each 
other. It is somewhat surprising that the squared Meissner mass is 1/3 of the squared Debye 
mass, as it is for the U{1)b case. Since the ratio of these two masses is the velocity of the 
Goldstone bosons, we flnd that the velocities of all Goldstone modes in our theory is equal 
to the speed of sound 1/V3 to the leading order of perturbation theory.^ 

4 Meson masses 

We now turn on flnite small quark masses and compute the resulting masses of the Gold- 
stone bosons. By introducing bare quark masses into the microscopic theory, we break the 



■^Similar results, v = vp/y/S, have been found by Bogolyubov and Anderson for phase waves in BCS 
superconductors and by Leggett for spin waves in superfluid ''He |18|] . 
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SU{3)l X SU{3)ji symmetry. This means that more terms are allowed in Ccs- However, 
constraints on the form of these terms can still be imposed if we note that the bare quark 
mass term 

AC = ^[M^R + h.c., (34) 

with M being the 3x3 mass matrix, could be made invariant under the SU{3)l x S'f/(3)/j 
if the matrix M is not passive under this symmetry but also transforms together with ipL,R 
in the following way, 

V'L^t/i^L, i^R^Un^PR and M -^ UlMUI (35) 

Any term in the effective Lagrangian, written as a function of S and M, must respect the 
extended symmetry (p!0|), (pS]). At large /i the U{1)a symmetry is effectively restored, which 
imposes an additional constraint on the possible form of mass terms in the chiral Lagrangian. 
Under the U{1)a transformation, the microscopic degrees of freedom transform as 

^L -^ e"^^L, ^R -^ e-'^-ijR (36) 



where C, is an arbitrary pure phase. Eq. ( p6D implies the following transformation law of X, 

Y and S, 

X^e-^'^X, Y^e^'^Y and E ^ e~^^^E, (37) 

The quark mass term (0) is invariant under the U{1)a symmetry if one requires that M 
transforms as 

M -^ e^'^M. (38) 

under the U{1)a- Thus, any mass term in the effective Lagrangian must be invariant under 
both the SU{?>)l x SU{3)r and the U{1)a symmetry extended by the transformation of M.0 
As the bare quark masses are assumed small, we want to construct the mass terms 
of lowest possible order in M. To the first order, one candidate is the mass term of the 
chiral Lagrangian at zero chemical potential, TrM'I'E. This term is invariant under the 
SU{3)l X SU{3)r, however it is not invariant under the U{1)a symmetry. Therefore, we 
must consider terms of higher order in M.0 

"^A similar symmetry argument allows one to fix the form of the effective Lagrangian in 2-color QCD to 
lowest order in fj. and ruq pAj . 

^The fact that A£cff — 0{AP) means that in the chiral limit M -^ the chiral condensate ('0'0) = 0, 
while {{ipip)"^) 7^ 0. Such a pattern of spontaneous chiral symmetry breaking has been discussed in [|l9[ . 
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Figure 3: The leading order diagram contributing to the shift of the vacuum energy by a 
small quark mass M. 

In the order M ^ we find that the following term is allowed by the symmetry^ 

A/:cff = -c-detM-Tr(M-^S)+h.c. (39) 

Another term allowed by the symmetry i^ 

A4g = -c' \li{MT)f + iTiMT^f] det S + h.c. (40) 

The coefficients c and c' of the mass terms from (|39|) and ( ^OD can be calculated by 
matching the shift of the vacuum energy as a function of S they induce in the effective 
theory: A£^vac(S) = A£eff + ^'^eff) to the perturbative calculation in the microscopic theory. 



The nontrivial, S-dependent, shift in the vacuum energy due to (|3^) is given to the lowest 

order by the diagram in Fig.^, which is similar to the second diagram in Fig. 0. Another 

diagram, similar to the first diagram in Fig |l|, also contributes to the shift of the vacuum 

energy, but its contribution does not depend on S.0 As a result, terms such as Tr(MS^M^S) 

and (TrME''')(TrM"^S), in principle, allowed by the symmetry, are absent to the order we 

are working. 

The calculation of the diagram in FigJ^ gives: 

51 + 321n2u2 ^ , 15-161n2u2 

c = -Ht and c = -Hr- (41) 

108 27r2 216 27r2 ^ ' 



^ One way of looking at this term is to realize that under S'C/(3)l both M and S transform as fundamental 
3-plets. We require two powers of M and one of S in order to satisfy the U{1)a neutrality. We can construct 
an SU{3)l singlet out of a product of M, M and E if we antisymmetrize with respect to the first index of 
each of these matrices (the index on which 5'f/(3)L acts). Antisymmetrizing also with respect to the second 
index we obtain an SU{3)l x SU{3)b. singlet: ea' b' d ^abcMaa' Mbb''Scc' , which coincides with ( p9| ) up to a 
constant. 

^ The alternative linear combination of the two terms in ( fJOf ) is equivalent to the term (|39|): 
[Tr(MSt)2 _ (TrM St)2] det E = 2 det M ■ Tr(M-iS). 

®This is related to the fact that such a diagram does not depend on the values of the gaps A on its two 
fermion lines, even if they are different. 
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The coefficient d is numerically very small: d jc ~ 0.027, and is related to the fact that 
one of the 9 gaps in the fermion propagator matrix (|23| ) is different from the other 8. This 



completely determines the mass term in the effective Lagrangian given by (^) and (|iO|). 
To write the (mass)^ matrix for the Goldstone bosons we expand the field S: 

"^'^H— h^"— -^i^^-' '''' 



2 



where a = 1, . . . 8, 9, A^ = \/2/3 and tt^ = vl j^^j f^i. Since the kinetic term is conventionally 
normalized we can read off the (mass)^ matrix, M\^ M.' ^^, from the mass terms (|39D, ( ^OD 
(except for the trivial rescaling of the r]' field): 

-A^^^vrV = -detM-Tr(M^U'^A^)7rV; 

^A^'^-^TT^ = C"{7rV[(TrAfAJ(TrMA5) + (TrM)(TrMA„\) 

+ (TrMA„MA6) + (TrM^A^Afe)] + 'ii'o'f [(TrM)' + (TrM^)" 

-2V6 7r'^7r^ [(TrAf)(TrMAa) + (TrM^Aa)] } ; (43) 

where, using (|32|) and (^T]), 

_ 2c 151 + 321n2 , _, 2c' 115-161n2 , , 

C = — = — ^ 1.578 and C' = — = ^ 0.04216. 44 

/2 3 21-81n2 /2 6 21-81n2 ^ ' 

Since these are pure numbers, the meson masses depend only on the bare quark masses, in 
contrast to the QCD vacuum where they also depend on the value of the chiral condensate. 

In order to understand the pattern of masses let us begin with the first term in (|43|), 
since the coefficient of this term, C, is much greater than that of the second term, C . 

With the quark mass matrix M = diag(m„, rud, ms), the 9x9 meson mass matrix ( ^3|) 
decomposes into a diagonal 6x6 matrix and a non-diagonal 3x3 matrix. The former gives 
rise to the mass formulas for vr^, K^, K^ and K^: 

ml± = C(m„ + md)'ms + . . . ; mj^± = C{mu + ms)md + . . . ; 

m|-o = C{md + ms)mu + • • • , (45) 

where ellipses denote contributions proportional to C . The remaining 3x3 matrix cor- 
responding to the 7r°, 7] and rj' mesons is not diagonal. The mixing pattern is easy to 
understand if we neglect the small difference between /,r and /^/. Then the 3x3 mass 
matrix has the simplest form in the basis mm, dd, ss, related to 7r°, t], t]' by the quark-model 
relations tc^ = {uu — dd)/^/2, r] = {uu + dd — 2ss)/\^, and rj' = [uu + dd + ss)/v^. In this 
basis \uu = diag(l,0,0), Aj^ = diag(0, 1,0), Ajs = diag(0,0, 1), and the 3x3 mass matrix 
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becomes diagonal. The mixing between Tr*^, 77, and 77', thus, is ideal, and the masses of the 
mixed states are 

m|„ = 2Cmdms + . . . , m|^ = 2Cm„m^ + . . . , m?^ = 2Cm„md + . . . . (46) 

As we shall see, however, the second mass term makes a significant correction to this ideal 
mixing, despite the fact that C ^ C, when m^ ^ rnu,d- 

The reader will notice that the mass ordering induced by the first term in (^) is com- 
pletely reversed compared to what one sees in the QCD vacuum. While this is quite surpris- 
ing, and seems unnatural under the hypothesis of continuity between quark and hadronic 
matter, it can be easily explained. The key point is that the mesons, which are the fiuctua- 
tions of S in Eq. (^), should be thought of as bound states of a triplet antidiquark X and an 
antitriplet diquark Y\ i.e. as qqqq states, rather than qq states. Consider, for example, the 
ss state. Now the s quark should be replaced by the ud antidiquark, and the s antiquark 
should be replaced by the ud diquark, so this meson is represented as udud. Since such a 
meson does not contain the strange quark, it is not surprising that its mass does not depend 
on nts. For all other mesons, one can write the quark structure as well, by making the 
replacement u —>■ ds, d —>■ us, s —^ ud. Since one replaces the heaviest quark by the lightest 
antidiquark and vice versa, the inverse mass ordering can be expected.0 It is important to 
note, however, that the mesons in the CFL phase have the same quantum numbers (up to 
mixing) as the mesons in the QCD vacuum. 

Although the coefficient of the second term in (|43|) , C, is almost 40 times smaller than 
C, it contributes significantly to the masses of some of the mesons, because the ratio of mg 
to rrici or rriu is also very large. The masses of the n and K mesons are given by (^. In the 
limit rus ^ m^^rf they reduce to: 

m1± ^ C {rriu + rriiijms; m|^o ~ (Cm„ + 4C"ms)ms; m|^± ^ (Cmrf + 4C"ms)ms;(47) 

where we treated both C/C and ms/md as large numbers of the same order of magnitude. 
We see that in the kaon masses the smallness of C jC is compensated by the greatness of 
nis/mu^d- The ordering of charged and neutral kaons remains inverted, however. 

The remaining 3x3 matrix for ttq, t] and r/' is rather complicated. To understand its 
properties it is helpful to consider the limiting case when m„ = irtd = 0. In this case 
the matrix simplifies, and one finds that two of the three states are massless (together with 
charged pions, in accordance with the Goldstone theorem and the breaking of SU{2) x U{1)a), 
while the third has a mass of order iris. The massless states are pure 7r° and (2v^vr^ — 7r^)/3, 
which is mostly 77'. 

^It has been suggested that a similar pattern arises for scalar mesons in QCD vacuum |2Ct|. 
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In order to get the feeling of the numbers involved, we substitute for the bare quark 
masses the values rriu = 4 MeV, rrid = 7 MeV and rris = 150 MeV. Performing numerical 
diagonalization of the mass matrix ( ^3]) we find the following spectrum: m^ = 117 MeV, 
m.,^0 = 53 MeV, m^/ = 30 MeV. The mixings of tt" with t] and rj' are on the level of a few 
percent, while the mixing between rj and 1]' is on the level of 20%. The numerical values for 
the masses (^ of pions and kaons are given by: m^± = 53 MeV, -m^o = 76 MeV, mx± = 72 
MeV. 

5 Conclusion 

In this paper we have shown that all parameters characterizing the dynamics of Goldstone 
bosons in the CFL phase, including the decay constants, masses, and the maximum velocities, 
can be reliably calculated in the weak-coupling regime of high densities. These parameters 
depend only on very few inputs — in fact, only the decay constants depend on the chemical 
potential, while the meson masses depend only on the bare quark masses and the velocities 
are given by a constant. This fact allows us to derive the full chiral Lagrangian without being 
dependent of the calculation of, e.g., the BCS gap, whose asymptotic behavior is known [0, 
but the exact numerical value of the prefactor is not [0, ||] . 

In contrast to the meson masses in the QCD vacuum, we found that, when the strange 
quark is much heavier than the u and d quarks, 5 mesons have masses of the order of m^, 
and 4 mesons have masses of the order {rriu^msY^'^. The lightest of those mesons is rj' with 
a mass around 30 MeV. The most surprising fact we found is the partially inverse mass 
ordering of the meson nonet. As we explained, this can be understood if we treat the mesons 
as qqqq states, rather than qq states. 

We must also discuss the region of validity of our treatment. The effective theory we 
are working with contains only meson modes, therefore, it is applicable only for energies 
smaller than 2A. In particular, a meson exists only as long as its mass is smaller than 2A, 
otherwise it would rapidly decay into a particle-hole pair. At asymptotically large chemical 
potential where A increases with /x, all mesons are stable with respect to this decay channel. 
At intermediate densities, the gap may actually drop substantially below 100 MeV in some 
range of the chemical potential 0. If the gap is small, some of the mesons may become 
unstable. For example, for our values of the quark masses, the rj meson disappears if A is 
smaller than about 60 MeV, and the lightest meson, r/', disappears only when the gap drops 
below 15 MeV. 

Another effect which becomes important when /i is not sufficiently large is the explicit 
breaking of the U{1)a symmetry by anomalous instanton-induced interactions. This effect 
gives a direct contribution (independent of the quark masses) to the mass of the rj' state, 
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and vn}!'^ contributions to the masses of the other mesons 0] . 

We can estimate the value of /i at which the transition from the normal QCD vacuum to 
the CFL phase occurs very roughly by comparing the value of /^r ~ 0.2/i in the CFL phase 
to its vacuum value /^r = 93 MeV. We find /itr ~ 500 MeV, which is reasonable and agrees 
with other estimates. 

It is straightforward to extend our calculations to the case of finite temperature. Let 
us discuss, qualitatively, the case of temperatures close to critical. The Meissner masses 
should decrease as the temperature increases and vanish at the critical temperature, since 
one expects the Meissner effect to be absent in the normal phase. On the other hand, Debye 
masses are nonzero in the normal phase, and so should not vanish at the critical temperature. 
Thus, near the phase transition, the velocities of the Goldstone modes become small, which 
is the manifestation of the critical slowing down of the relaxation of the condensate near the 
phase transition. The decay constants remain finite, while the masses go to zero as yuA/Tc 
near the phase transition. It would be interesting to study possible consequences of this 
behavior. The inverse mass ordering and anomalous lightness of mesons, in particular, of 



the 77' [21 1, in the CFL phase may also significantly affect their production should such phase 



be accessible in heavy ion collisions. 
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